Vacuum degeneracy of a circuit-QED system in the ultrastrong coupling regime 
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We investigate theoretically the quantum vacuum properties of a chain of iV superconducting 
Josephson atoms inductively coupled to a transmission line resonator. We derive the quantum 
field Hamiltonian for such circuit-QED system, showing that, due to the type and strength of the 
interaction, a quantum phase transition can occur with a twice degenerate quantum vacuum above a 
critical coupling. In the finite-size case, the degeneracy is lifted, with an energy splitting decreasing 
exponentially with increasing values of g 2 N 2 , where g is the dimensionless vacuum Rabi coupling 
per artificial atom. We determine analytically the ultrastrong coupling asymptotic expression of the 
two degenerate vacua for an arbitrary number of artificial atoms and of resonator modes. In the 
ultrastrong coupling regime the degeneracy is protected with respect to random fluctuations of the 
transition energies of the Josephson elements. 

PACS numbers: 



Circuit quantum electrodynamics (circuit-QED) is a 
very fascinating topic for fundamental condensed mat- 
ter physics, quantum optics and quantum information. 
In superconducting circuit-QED systems, it has been 
possible to implement on a chip the celebrated Jaynes- 
Cummings model by strongly coupling a superconducting 
artificial atom to a bosonic mode of a microwave trans- 
mission line resonator [HQ and to perform quantum log- 
ical operations with two qubitsjU. So far, experimental 
manipulation of quantum states in such circuit-QED sys- 
tems has dealt with excited states. In these systems, the 
quantum ground state is non-degenerate and no infor- 
mation can be stored or processed by using only the vac- 
uum of the circuit-QED system. In principle, by a judi- 
cious choice of their components, superconducting quan- 
tum circuits can give rise to Hamiltonians, which cannot 
be achieved in atomic cavity-QED systems: in particular 
it may be possible to tailor the relative amplitude and 
the form of the interaction terms for the generation of 
interesting and controllable quantum vacuum properties. 

Here, we present a rigorous quantum field derivation 
showing that it is possible to obtain a vacuum degeneracy 
of a circuit-QED system by using a chain of Josephson 
junction atoms inductively coupled to a transmission line 
resonator. A quantum critical coupling occurs in such a 
circuit-QED system thanks to both the type and ultra- 
strong size of the interaction obtainable with the induc- 
tive coupling scheme. In the case of a finite number N of 
artificial atoms, a degeneracy lifting occurs, with an en- 
ergy splitting dramatically decreasing as cxp(— g 2 (3(N)), 
where g is the dimensionless vacuum Rabi coupling per 
atom (i.e., vacuum Rabi frequency divided by the tran- 
sition frequency wj?) and /3(A) depends quadratically on 
A. We present the asymptotic formula of the two degen- 
erate vacua in the ultrastrong coupling limit . Moreover, 
we show that the degeneracy is protected with respect 
to random site-dependent fluctuations of the Josephson 
transition energy. 
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FIG. 1: A chain of N Josephson atoms ("F" stands for 
fluxonium [Ij]) are inductively coupled to a transmission line 
resonator. By tuning the external magnetic flux, the flux- 
dependent potential for each fluxonium has a symmetric dou- 
ble well structure with two states |0) and |1) (with energy dif- 
ference TvjOf) forming the two- level system (parameters used 
for the inset: Ej/E 0j = 3, Ej/E Lj = 20). 



A sketch of the proposed system is depicted in Fig. [TJ 
namely a chain of A identical artificial two-level atoms 
in a transmission line resonator. Each artificial atom 
('fluxonium' is made of a Josephson junction cou- 
pled to inductances and an external magnetic flux. Here, 
we consider a scheme where each fluxonium is inductively 
coupled to the resonator. The fluxonium is known be free 
from charge offsets moreover, the inductive coupling 
can produce extremely large coupling even with a single 
artificial atomQ. In the case of a chain, the Hamilto- 
nian is H = ^ . Hj where each cell of size a is labeled by 
the index j and is located at the position Xj. One can 
effectively model the resonator as a sequence of indue- 
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tances L r = al r and capacitances C r — ac r [6[, where l r 
(c r ) is the inductance (capacitance) per unit length. The 
Hamiltonian of each cell reads: 
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2L X 



%^ + iE Cj (N^~E J cos(^ + ^ ext )) ,(1) 



where N and <f> are the number and flux operators for the 
resonator elements and Josephson junctions (V stands 
for resonator; 'J' for Josephson junction). The charging 
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energies are Eq t = -%q- and Eq, = jfr" . By applying 
Kirchoff's laws and by taking $ ea;t = 7TjS we find H = 



Hres ~t~ E[p -\- H COU pH n g 1 
N 



with: 



Hres = £ ^C r (Nr j ) 2 + E Lr ^ ~f 



3 = 1 



iv , „j >2 

Hp = ^EcAN'jf + E Lj ^f- + Ej cos(#) , 
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Hcouphng = G( M - £r ( 2 ) 

3=1 

where we have introduced the dimensionlcss fluxes tp^ = 
-^44, = x^J ana - tne inductance energy con- 



stants arc E^ r 
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The magnitude of the cou- 
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(JL) 2 ( L t +L r \ 

\2e' [ -L 1 L r + L 1 L 2 +L 2 L,J 

pling constant is G = 

Hamiltonian H res describes the transmission line res 
onator with a renormalized inductance per unit of length 



-L 2 + - 



" r ~ " r L t +L 2 

ductances in each fluxonium 
ment in Ref. 



accounting for the additional in- 
By following the treat- 
the position-dependent flux field is 

fa) = iT, k >i £y/Wf k ( x ) ( &k ~ a ^ wherc a fc is thc 
bosonic creation operator of a photon mode with en- 
ergy hjj k = ^f-^J%E Cr E Lr . The spatial profile of the 
fc-th mode is fk(x) = — \/2/dsin(^p) for k odd, while 
fk(x) - 



2/dcos(-^p) for k even, d being the length of 
the one-dimensional resonator (in the following, we will 
consider d = Na). The site-dependent fluxes are simply 
given by cj>> r = (f>(xj). 

Thc Hamiltonian Hp describes the sum of the ener- 
gies of the fluxonium atoms. By properly tuning the 
external magnetic flux, it it possible to obtain a symmet- 
ric flux-dependent potential energy, as shown in Fig. [TJ 
with a double well structure. Due to the strong anhar- 
monicity of its energy spectrum, the fluxonium can be 
approximated as a two- level system, when Ej Elj- 
We call the two first eigenstates of the j-th fluxonium 
and \l)j and we introduce the raising operator 



as |0)j 



11X01,- and &_ d 



\0)(l\j. By using thc 



Pauli matrix notation, we have a x 
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and 



cfy.j = i(o-+,j - and a z 



2 °+,3 a +,j 



1. Leaving 



aside a constant term, we then have Hp = . Tuvp^azj, 
where fuup is the energy difference between the two states 
|0) and |1). By considering only the two-level subspace, 
the Josephson junction flux has the form 
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where ipoi — ^ for typical parameters ( sec Fig. 1). 

As it will be clear in the following, it is conve- 
nient to introduce excitation creation operators St = 

N^f=i^fk( x j)v+,3 for 1 < fc < iV — 1, where 

Af k {xj) = cos( fc " ( ~y 1+j) a) for k odd , and Af k (x 3 ) = 



sin( 



fc7r(- 



-+3). 



for k even. Note that the collective 



operator b\ is a linear superposition of the excitation op- 
erators in each fluxonium with an amplitude depending 
on the spatial profile of the flux field of the resonator. In 
order to get a unitary transformation, it is also necessary 
to introduce the operator V N = -== ^V^+J- ^ n 

the following, we will consider only the resonator modes 
1 < k < N, because, in the conditions we are consid- 
ering, the higher order (Bragg) modes are energetically 
well off-resonant. Hence, we get the following effective 
Hamiltonian: 



H 



Kk<N 



(4) 



where = (ah, bk, at, b\) T with thc Bogoliubov di- 
agonal metric rj = diag[l, 1,-1,-1], and the matrix: 
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The coupling between the annihilation operators d k , b k 
and the creation operators are due to the antiresonant 
(non-rotating wave) terms present in the coupling Hamil- 
tonian H coupling- The collective vacuum Rabi frequency 
reads for 1< k < N - 1 
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hYl k = G—tpoi sm — — 
n za ui k 
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( and for k = N, m N = Gf^^f). 
Notice that each fc-mode of the resonator is coupled 
only to the collective matter mode with the same spa- 
tial symmetry and H = X^fc^fc- Hence, the eigenstates 
are products of the eigenstates corresponding to the k- 
subspaccs. Thc effective Hamiltonian in Eq. (j4|) has 
been obtained by assuming that the operators b\. are 
bosonic, i.e. [&&,&£] ~ 1, an approximation working in 
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the limit N 1. The excitation spectrum of the col- 
lective bosonic modes depends on the eigenvalues of the 
matrix M. k . A crucial property is given by the deter- 
minant Det(AAk) = WkU!p(uikUiF — 4f2?), which vanishes 
when the vacuum Rabi frequency equals the critical value 
Q£ = ^"^'" F | implying that two of the 4 eigenvalues of 
Aik are exactly zero. For Qj. > fi*, two of the 4 eigenval- 
ues of the matrix M. j. becomes imaginary, manifesting an 
instability of the normal, non-degenerate, quantum vac- 
uum phase. 

This is reminiscent of quantum phase transitions with 
Dicke-like HamiltoniansQ, where at the quantum criti- 
cal point there is a gapless bosonic excitation. Note that 
Dicke-like Hamiltonians arc usually obtained by dropping 
the so-called A 2 -term, that is the term associated to the 
squared electromagnetic vector potential term (To| . 
However, because of the magnitude of the A 2 -term, a 
system with (ultra)strong light-matter coupling does not 
necessarily have a quantum critical point For exam- 
ple, in the celebrated Hopfield modelfTJj for dielectrics a 
quantum critical point does not exist even if the coupling 
can be ultrastrong 13lll4|. In our case, the quantum field 
Hamiltonian in Eq. ([2]) is the complete Hamiltonian for 
the proposed superconducting system and no term has 
been omitted. In particular, the analogous of the A 2 
term in the present system is given by extra terms pro- 



portional to ((fil 



_1 ) 2 obtained from Eq. (1) after 



substitution of the expression for the flux ij>> x obtained by 
Kirchoff's laws [Til]. Here, these terms are fully included 
and contribute to the expression for the resonator renor- 
malized inducting energy Ez, r , which does depend on L\ 
and L2. 



For the case of finite number of fluxonium atoms ./V and 
finite number of modes N m , we have performed numeri- 
cal diagonalizations of the circuit-QED Hamiltonian. As 
shown in Fig. 2, for increasing coupling the energy of 
the first excited state converges towards the energy of 
the ground state, hence a twice degenerate vacuum is 
obtained in the ultrastrong coupling limit. As shown in 
the inset, for a given value of N, the energy splitting ex- 
ponentially decreases with increasing vacuum Rabi cou- 
pling. The finite-size scaling properties are shown in Fig. 
3(a), where the energy splitting is plotted as a function of 
q 2 , where q = 3°=' — is the dimensionless vacum Rabi 

, VNuj k = 1 

coupling frequency per fluxonium. Our numerical re- 
sults show that the energy splitting scales exp(— g 2 (3(N)), 
where f3(N) w 2N 2 (see inset of Fig. 3a). Hence, a per- 
fect degeneracy is obtained either in the thermodynam- 
ical limit (N — > +oo) or for g 3> 1. As shown later, 
it is possible to have g > 1 in realistic superconduct- 
ing systems, hence a negligible splitting can be achieved 
with a relatively small N. In the ultrastrong coupling 
limit ( °J;= 1 — > oo), we have derived [li} an analytical ex- 
pression for the two degenerate ground states by taking 




FIG. 2: First 30 energy eigenvalues versus dimensionless vac- 
uum Rabi coupling for TV = 5 fluxonium atoms, N m = 3 
resonator modes (cjf — Wfc=i)- Upper inset: the difference 
between the energy eigenvalues and the ground state energy 
is plotted. Lower inset: normalized energy difference (log 
scale) between the first 2 quasi-degenerate levels versus the 
dimensionless coupling. In the ultrastrong coupling limit, the 
two ground states are excellently approximated by the analyt- 
ical formula in Eq. J7J (98% overlap for the largest coupling 
here considered). 



into account an arbitrary number N m of modes for the 
resonator. It is convenient to introduce the x-polarized 
states (eigenstates of cf x> j), namely \+)j = + 

and = ^7|(|l)i — We have found that in the 

ultrastrong coupling limit (where Hp is dominated by 
H res + H coupling) the asymptotic expression for the two 
degenerate vacua |G+) and |G_) is: 



\G±)=C G n j \±) j ®IL ko e 



±<7 



gV2 



~ al) \o) ko 



n fe jo)i 



(7) 



with Cq a normalisation constant, k a (k e ) standing for 
the odd (even) k values for the resonator modes. Eq. 
([7]) shows that the two degenerate ground states are the 
product of a 'ferromagnetic' state for the chain of arti- 
ficial atoms times coherent states for the different res- 
onator modes. Importantly, the two orthogonal ground 
states have opposite polarization of the pseudospins and 
opposite phases for the coherent states. Due to the mode 
spatial symmetry, in \G±) the even k e resonator modes 
are empty. The analytical expression for the two vacua 
excellently agrees with the numerical results. 

It is interesting to see how the degeneracy is affected by 
the presence of an additional Hamiltonian term H pert = 
Y^j ^hAjCTzj, describing a site-dependent random fluc- 
tuation of the fluxonium energies. Interestingly, we have 
found numerically (see Fig. 3b) and analytically 
that the average splitting < S > and its standard de- 
viation a = \/ < 5 2 > — < 5 > 2 have the same expo- 
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limit of adiabatic changes of the Hamiltonian, thanks to 
the degeneracy, it may be possible to create non-abelian 
Berry phases and control quantum superpositions in the 
ground state subspace (a sort of vacuum qubit). The 
system studied here appears promising for the observa- 
tion of quantum phase transitions and the manipulation 
of quantum vacua in circuit-QED. 

We are grateful to M.H. Devoret for many discussions 
and for suggesting us to consider the nuxonium atom. 
We thank B. Dougot, I. Carusotto, S. De Libcrato for 
discussions and M. Bamba for numerical advice. 



FIG. 3: Normalized energy splitting (log scale) between the 
two quasi-degenerate vacua versus g 2 for different values of N 
(2, 3, 4,.. .11), where g is the dimensionless vacuum Rabi cou- 
pling frequency per artificial atom, (a) Results for the case 
of identical nuxonium atoms with same transition frequency 
ujf- 5/ujf decreases as exp(— f}{N)g 2 ). Inset: /3(N) versus 
N 2 . (b) Averaged splitting < 5 > in presence of a random 
distribution: ujj t F = <^f + Aj = Uf(1 + 0.5£j), where j labels 
the site and is a gaussian variable with variance equal to 
1. The results have been averaged over 100 disorder configu- 

v< s 2 > 



rations. The standard deviation a 
the same exponential dependence (see inset) 



< S > 2 has 



nential dependence as the splitting in the ideal case of 
identical Josephson elements. By having N and/or g 
large enough, the effect of a disorder of given ampli- 
tude can be made arbitrary small. This occurs because 



(G±\H™ rt \G±) 



(G±\H™ t \G T ) = with m<N-l, 



i.e., such a perturbation is zero up to the iV-th order per- 
turbation theory, leading to a protected degeneracy 17|. 



The present system can indeed have a very large con- 
trollable coupling. For ujp = Wfc=i, we have 
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SUPPLEMENTARY ONLINE MATERIAL FOR THE PAPER : VACUUM DEGENERACY OF A 
CIRCUIT-QED SYSTEM IN THE ULTRASTRONG COUPLING REGIME 

In Section A of this Supplementary Material, we give a detailed proof of our analytical expression for the two 
degenerate vacua in the ultrastrong coupling regime. In Section B, we show the analytical calculation of the degeneracy 
splitting in the ultrastrong coupling limit. We show also useful details about the degeneracy with respect to site- 
dependent fluctuations of the Joscphson atom transition frequencies (energy disorder). 



A. ULTRA-STRONG COUPLING LIMIT: ANALYTICAL APPROACH 

Let us consider the ultrastrong coupling limit ( -^y 1 — > +00) , by taking into account an arbitrary number of 
modes N m for the transmission line resonator and a finite number TV of fluxoniums qubits. The following treatment 
is general and can be done for any kind of mode spatial profile Af k (x). In the ultra-strong coupling regime, the 
bare Hamiltonian Hp of the pseudospins (the Josephson artificial atoms) can be treated as a perturbation of the 
Hamiltonian H res + H coup ii ng , where H res is the Hamiltonian of the bare transmission line and H coup u ng is the 
interaction part between the resonator and the artificial atoms. As written in the letter, it is convenient to consider 
the basis (|— ), \ +))j, the two eigenvectors of the operator & x ,j- The Hamiltonian can be rewritten in the form: 

N 

Hres + H coup li n g = E ^4°* + E E ^ fc V A ^ k ^ ) ^ k ~ 4 ) ( ( I +) (+ \j ~ I ~) (~ 1 3 ) ( 9 ) 
k=l..N m k=l..N m j=l 

Let us introduce the subspace J r s < generated by the states II j 1 ) ® l^res) where the Q G { — , +} denote a given 
pseudospin configuration sequence for the chain, while |\& res ) describes a generic state of the resonator bosonic 
quantum field. Since there are TV two-level systems, then we have 2 N subspaces ^F{\^ u =1 N , each for every pseudospin 
configuration. 

It is clear that by applying H res + H coup u ng on a state \ip) <S -^V I Cj > > =1 nj ^he result is a state belonging to the same 
subspace (i.e., H res + H coup u ng conserves a given pseudospin configuration). Now, if we consider a given pseudospin 

configuration sequency for the chain, we can define -0f c = J2j \fj^^fk( x j)ns < (j) where iis c (j) = 1 if l^o) = 
and where ns c {j) = — 1 if \Sq) = we have that the Hamiltonian on the (J^Sc) can be written as: 

Hris + H S j uplmg = E (tuAa k +ihn k (a k -al)i> S k <)= £ M4 + - ^) - Av£ e ) a ). 

(10) 

So if we introduce the shifted boson operator a k < = a k — i^-tp k c , we have the same bosonic commutation relations 
between a fc c and (ct fe c )^ and the Hamiltonian terms read 

h s <= E ^(af^a? - E ^k f (n) 

fe=l..JV m k=l..Nm k 



So, on this subspace, the fundamental state \Gs ( ) has the energy Eq s = — J2k=i N m ^■T^(' l l J k < ) ■ Moreover, we 
have: 



a S k ( \G S( )=0 Vfc<TV„ 



This relation implies that: 



\G S< ) = n^)®!!^—^—^ 1 ^ a - ] \0) k (12) 



In order to find the ground state for the complete Hilbcrt space, we have to determine the pseudospin configuration 
minimizing - £fe=i..iv m 
called Q the quadratic form 

K £lk=i k 



minimizing - Efe=i..jv m ^Sr(V J f C ) 2 = -2.9 2 (fiwfc=i) Y,j,j> Vs t ij)Q(j, j')vs t {j') where g = J^J^ =i and where we have 

Q(j,j')= E A/ fe (x,)[(-^) 2 i]A/ fe ^.). 



k=l..N n 



6 



For a given number of modes N m and a given spatial profile Afk(x) and position (xj)j—i..N, we derive the form and 
we find beyond the 2 N configurations S^, the one which minimizes — Ylj,j< ^s c {j)Q(j, j')ns ( {j')- In fact the double 
degeneracy of the spectrum appears also here because at any configuration corresponds, an opposite one S$> (for 
which ns^U) — ~f jt S((j) V?) with same energy. With our particular profile , the configurations of minimal energy 
are the two ferro-magnetic ones: fis ( (j) = +1 Vj and fJ-s c ,(j) — — 1 Vj. So, to conclude, the two fundamental states 
we derived are : 



\G±) = n,|±) ® n fc e" '" fc t fc ' e (^v- t s± 4)| 0)fc = CaU^j ® n feo e ±(fc ^' ; <^» Qlj |o) feo ® n fc jo) fee (13) 

where Gg = n^e 2 and with k Q (k e ) standing for the odd k values for the resonator modes. 



B: DEGENERACY SPLITTING AND PROTECTION WITH RESPECT TO SITE-DEPENDENT 
ENERGY DISORDER IN THE ULTRASTRONG COUPLING LIMIT 

Now let us consider the effect of the bare Hamiltonian for the artificial atoms, namely Hp = ^F^jLi ~ 

hjjp Y^j=i \j + I - where |±}j are the eigenstates of a x j. In the ultrastrong coupling limit, Hp acts as 

a perturbation of H res + H coup u ng . In the finite-size case, this produces a degeneracy splitting for the two vacua |G+) 
and \G-). In presence of TV artificial atoms, the effect of Hp is zero up to the N-th order in perturbation theory. 
As shown by the numerical results reported in our letter, the splitting decreases as ~ q-p\ n )9 w ith f3(N) ~ 2N 2 
(see inset of Fig. 3a). Here we show in detail for the case N = 2 that this exponential dependence can be found 
analytically. Moreover, the same result occurs in presence of a random site-dependent fluctuation of the artificial 
atom transitions energies (energy disorder). We finally discuss the N> 2 case. 



Degeneracy splitting for N=2 fluxoniums in the ultrastrong coupling limit 



Let us consider the term tuj F h(\+)(— |i + |— that is the bare Hamiltonian term associated to the first artificial 
atom. Taken alone, such a term does not lift he degeneracy since it does not couple directly |G+) and |G_) at any 
order. However, it produces a mixing of the vacuum states |G+) and |G_) with the excited states of H res + H coup n ng . 
Accordingly, we can introduce the following states 



\G+) 
|G_) 



IG, 



(n,- + |^|G H 
Eg+ — E n — [ 



n,- 



. 1. (n,+ -|of|G_). 



(14) 
(15) 



where n = (ni,ri2, Hjv m ) and the states |n, — h) 



((H=iV) ni ---((H=Njr Nm \G-+) stand for the 



excited eigenstates of H re 



En,- 

J 



y/m\-n Nm ] - vv " , fc= 1/ ' ■■■^^k=N„ 

HcoupUng with pseudospin configuration { — +}. Their corresponding eigenenergies are 

7iT+\\ = 



Hk=i n n kf>^k + Eg_ + ■ They are obtained by applying the shifted photonic creation operators (a k ) 



■ i—tp k + for the mode numbers k = 1 to k = N m . Anagolous definition holds for |n, H — ), corresponding to the 



pseudospin configuration {H — } . 

Now, we can take care of hiop^a^ = fkOp^a^\+)(— [2 + | — X+I2), that is the contribution of the second pseudospin. 
The energy splitting between the vacua |G+) and |G_) at the N = 2 order perturbation theory is 



5 = uj F {G-\cr z 2 \G + ) 



ojp f ^ nu P {G +\ai\n, +-)(n, + - \oj\G-) 



Eg, 



E n _j — 



E 



^<G_|af|n,-+}<n,- + |^|G H 



En, — h 



E 



E G - 

hu F {G + \a%\n, +-)(n, 



E G . 



Wf\G-) 



E n _j — 



.(16) 



We can now consider the effect of only one photonic modc( N m = 1). In that case, the result is simplified because 
ipk=i = ~ Afk=i{xi) + Afk=i(x2) = 0, hence the excited states \n, H — ) are simply the unshifted Fock states \n). 



7 



Then, the energy splitting reads: 



Wfc=l 



- 2 9 ia *= 1 )|n)(n|e- 2flia *=i)|0) 



4g 2 + n 



Wfc=l 



^L P -4s 2 V 



(-% 2 ) r 



( V + n)n\ 2uj k=1 




(17) 



In the last derivation we used the identity V'jta = A/fc = i(a;i) + Afk=i(x2) = \/2 = — 4>k=x an( ^ Eq + = Eq_ = 
— 4huik=ig 2 ■ Note that we have verified that this analytical expression is an excellent approximation of the exact 
numerical results (for the parameters used in the figures of our manuscript the analytical approximation differs less 
than 10% from the numerical results). 

Degeneracy splitting for N — 2 in presence of site-dependent energy disorder in the ultrastrong coupling limit 

Now, wc wish to consider a configuration where the artificial atoms have not the same energy. In Fig. 3b of our 
manuscript, we have reported numerical results showing that a disorder- induced degeneracy splitting dramatically 
decreases with the coupling per fluxonium g and size N. In particular, we have shown that there is the same 
exponential dependence as in the absence of disorder. This can be proved analytically for the N — 2 case. We 
consider 

Hp = ^ t'-l.j-f:.! ( 18 ) 

J'=l,2 



Following the same steps as in the previous derivation wc find: 



U}p 2 



E 



7KJF,i(G + |o-f|n, H — )(n, + - |of|G_) u F:2 u F ,i 



E n j — 



2 LU k=1 



(19) 



If we now consider different disorder realizations and we average over the configurations, we get an average splitting: 



< 8 >c 



2wfc=l 




-8g 2 



(20) 



Taking u)p^ = cop + Aj with A^ a random variable with zero average (< Aj >= 0), variance < A| >= A 2 and such 
that < AiAj^i >= Vi, j , then the averaged degeneracy splitting < <5 > and the standard deviation a read: 

2 

(21) 



\J < 5 2 > - < 8 > 2 




< 8 > 



(22) 



Hence, the averaged splitting is equal to the splitting without disorder for the case of identical Joscphson elements. 
The standard deviation of the splitting depends on the normalized disorder amplitude A/cop, but it has the same 
exponential dependence as < 8 >. Hence, in the ultrastrong coupling limit (g> 1), the effect of disorder can be made 
arbitrarily small. 

Note that a similar protection occurs with respect to local noise sources proportional to o y .j. 



Degeneracy splitting for N > 2 fluxoniums 



We now consider the general case N > 2 and Josephson artificial atoms with site-dependent energy, i.e. Hp = 
Ylj=i n ^F,j^z,j- In this general case, Hp couples the two degenerate vacua only at the iV-th order in perturbation 
theory. For finite values of N and <?, the energy splitting is given by the following expression: 



n 



j=l..N 



-) 



E E 



hS~ 2 

o-GSjv ni,n 2 ,...n N _i 

( G +k£(l)l n l;<r(l)(S + ) )( n l><T(l)(S + ) l cr a(2)l n 2) ( r(2)(<T(l)(S + )) > •••• ,<t(JV-1)(...<t(1)(S + )) \ a a(N)\ G ~) , 



(E G+ 



^ n lu(l)(S + ) )(^H.i7(l)(S + ) ^n2, I (2)(„(i)( S+ )))-(^nN- 



2,„(N-2)(.. CT (1)(S + )) 



ni>J-l>o-(w-l)(. 



(23) 
(24) 
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where is the set of permutations of {1..N}, S+ is the {+ + ...+} configuration, and where <r(m)(...a(l)(S+)) 
stands for the pseudo-spin configuration in which the a(l) th , a(2) th ... a(m) th pseudo-spins have switched from + to 
-. The last expression contains all the excited states of every pseudo-spin configurations, with their energies at the 
denominator. In fact, in the ultrastrong coupling limit, the denominator will give a polynomial contribution to the 
splitting proportional to (fojJk=i) N ~ l ■ Hence, we get 



n 



J=l. .AT I 



E E 



5 ~ 2ojfe = i . 

crGbw ni,n 2 ,...nisi_i 

(G , + |cr^ (1) |ni, .( 1 ) (5+ ) )(ni, cr ( 1 )( S+ ) |cr^( 2 )|n 2 , (T (2)(o-(l)(S + )) >----( n N-l,<j(Ar-l)(...cr(l)(S + )) Wa(N) \@-) 



2w k=1 N\ 



n 



3=1. -N 



LOp. 



2w k 



3 -) 



(G+|n i=1 ..^< (i) |G_) = 2w fe=1 iV! 

= 2u k=1 N\ 



n 



3=1. -N 



2w k= 



(G+|n i=1 ..^^|G_) 



Mj=i..N{—. > 



The averaged degeneracy splitting < S > and the standard deviation a over the disorder configurations read: 

. , -4g 2 ^ 1 

< o >~ 2u) k= iN<( ) e ( 2n' - m *e 

2wfc = i 



a = y/< S 2 > - < 5 > 2 ~ 2 Wfc=1 7V!(-^-) N e n; W ^ 

2w/c = i 



A 



(1 + ( — ) 2 ) N - 1 ~ \/iV— < (5 > 



So, at resonance, and keeping only the dominant term, 

log(!A^) = -/3(iV) 5 2 



-4g 
sin(^) 



V V 1 

JU 2 « 



i<fc c <Ar„ 



A-5 



sin (25r) : 



Si<fc <jv F" ^ 2.1iV 2 ViV > 2, which agrees with the numerical results . 



(25) 



(26) 



(27) 



(28) 



where 1.67V 2 < 

The standard deviation a of the degeneracy splitting does depend on the normalized disorder amplitude A /top, but 
it has the same exponential dependance than the average splitting 5. Hence, the effect of disorder can be made 
arbitrarily small when g and/or N are large enough. 



